Draft version November 19, 2009 

Preprint typeset using I4TpX style emulateapj v. 03/07/07 



UNDERSTANDING THE GEOMETRY OF ASTROPHYSICAL MAGNETIC FIELDS 



On 
O 
O 

> 
O 

OS 



< 

O 

6 



> 
os 
o 

Os 

m 



Os 
O 



Avery E. Broderick & Roger D. Blandford - 

1 Canadian Institute for Theoretical Astrophysics, 60 St. George St., Toronto, ON M5S 3H8, Canada; aeb@cita.utoronto.ca 
2 Kavli Institute for Particle Astrophysics and Cosmology, 2575 Sand Hill Rd., Menlo Park, CA 94309, USA 

Draft version November 19, 2009 

ABSTRACT 

Faraday rotation measurements have provided an invaluable technique with which to measure the properties 
of astrophysical magnetized plasmas. Unfortunately, typical observations provide information only about the 
density-weighted average of the magnetic field component parallel to the line of sight. As a result, the magnetic 
field geometry along the line of sight, and in many cases even the location of the rotating material, is poorly 
constrained. Frequently, interpretations of Faraday rotation observations are dependent upon underlying models of 
the magnetic field being probed (e.g., uniform, turbulent, equipartition). However, we show that at sufficiently low 
frequencies, specifically below roughly 13(RM/lradm" 2 ) 1/ ' 4 (B/lG) 1 / 2 MHz, the character of Faraday rotation 
changes, entering what we term the "super-adiabatic regime" in which the rotation measure is proportional to 
the integrated absolute value of the line-of-sight component of the field. As a consequence, comparing rotation 
measures at high frequencies with those in this new regime provides direct information about the geometry of 
the magnetic field along the line of sight. Furthermore, the frequency defining the transition to this new regime, 
fSA, depends directly upon the local electron density and magnetic field strength where the magnetic field is 
perpendicular to the line of sight, allowing the unambiguous distinction between Faraday rotation within and in 
front of the emission region. Typical values of j^sa range from 10kHz (below the ionospheric cutoff, but above 
the heliospheric cutoff) to 10 GHz, depending upon the details of the Faraday rotating environment. In particular, 
for resolved AGN, including the black holes at the center of the Milky Way (Sgr A*) and M81, vsa ranges from 
roughly 10MHz to 10GHz, and thus can be probed via existing and up-coming ground-based radio observatories. 
Subject headings: polarization - radiative transfer - radio continuum: general - magnetic fields - turbulence - 
Galaxy: center 



1. INTRODUCTION 

Magnetized plasmas are common in astrophysics, playing 
central roles in objects as diverse as galaxy clusters, the sites of 
star formation, the solar corona and accreting black holes and 
their ultra-relativistic outflows. Despite their importance, how- 
ever, there are few observational tools for assessing the physi- 
cal conditions within them. This is especially true for strongly 
ionized or diffuse regions, in which molecular line emission 
is absent. In these the best evidence for significant magnetic 
fields, beyond theoretical arguments, is due to Faraday rotation 
observations of intrinsically polarized sources. 

First discovered in optically active crystals (Faraday 1846), 
Faraday rotation has subsequently become one of the primary 
methods for measuring the strength and geometry of astrophys- 
ical magnetic fields. This has been primarily via the determina- 
tion of the rotation measure, defined in terms of the polarization 
angle, <&, by 

RM=-^ ~8.12x 10 5 / n.B-c/^radm" 2 , (1) 

where n e , B and t are measured in cm -3 , G and pc, respectively. 
Thus, RM's provides a line-of-sight averaged measurement of 
the line-of-sight magnetic field, weighted by the local plasma 
density. Measured values for the RM range from nearly zero to 
nearly 10 6 radm" 2 . 

Rotation measure observations have played a critical role 
in determining the magnetic field strengths in the centers of 
galaxy clusters (e.g., Ge & Owen 1993, 1994), Galactic mag- 
netic field strength and distribution (e.g., Men et al. 2008; 
Noutsos et al. 2008), magnetic field strength and structure 
within ultra-relativistic black hole jets (Zavala & Taylor 2004; 
Stirling et al. 2004; Miller- Jones et al. 2005; Brocksopp et al. 



2007; Kharb et al. 2009), magnetic field strengths near, and 
accretion rate of, the supermassive black hole at the center of 
the Milky Way and other nearby low-luminosity active galac- 
tic nuclei (AGN) (e.g., Agol 2000; Quataert & Gruzinov 2000; 
Brunthaler et al. 2006; Macquart et al. 2006; Marrone et al. 
2007), and even the solar corona magnetic field (e.g., Mancuso 
& Spangler 2000). 

However, there are a number of significant ambiguities in 
Faraday rotation studies. The first is the degeneracy between 
density and magnetic field strength (i.e., high densities and 
weak fields vs. low densities and strong fields). The second is 
the degeneracy between weak large-scale fields and strong tan- 
gled fields. In light of numerical simulations of magnetohydro- 
dynamic turbulence, and the observation of supersonic (though 
sub-Alfvenic) turbulence in Galactic molecular clouds, these 
make interpreting measured RM's highly dependent upon mod- 
els for the rotating medium. Since in many cases it is unclear 
if the Faraday rotation is occuring in situ or in a foreground 
region distant from the polarized source, there are large uncer- 
tainties in the our understanding of the physical conditions in 
the source. 

Here we show that at low frequencies Faraday rotation 
qualitatively changes, with the RM becoming proportional to 
f n e \B ■ d£\. Comparisons between RM's measured within this 
"super-adiabatic" regime and at high frequencies can probe the 
line-of-sight geometry of the magnetic field 1 . Furthermore, 
the frequency at which a given source transitions between the 
standard and super-adiabatic regimes depends upon the local 
plasma properties at magnetic field reversals (when the mag- 

1 Analogous phenomena include the "MSW" effect in neutrino astrophysics 
(Bahcall 1989) 
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netic field is orthogonal to the line of sight), and can thus dis- 
tinguish between different candidates for the site of the rotating 
medium. 

In Section 2 we review the plasma processes responsible 
for Faraday rotation and describe the super-adiabatic regime. 
Section 3 deals with the observational consequences unique to 
super-adiabatic Faraday rotation generally, while Sections 4 & 
5 concern specific classes of sources in which ordinary Faraday 
rotation has been observed. Our conclusions and implications 
are summarized in Section 6. The mathematical details of an 
ab initio determination of the radiative transfer regimes is pre- 
sented in the Appendix. 

2. FARADAY ROTATION & ADIABATIC PLASMA MODE 
PROPAGATION 

Faraday rotation is commonly presented in terms of the dif- 
ferent phase velocities of the circularly polarized eigenmodes 
of magnetized plasmas. Specifically, a linearly polarized, 
monochromatic, incident wave is decomposed into the two 
plasma modes, which subsequently propagate independently, 
accruing a net phase difference, A$ = J Akd£ (where Ak(uj) 
is the difference between the wave-vectors of the two modes), 
as a consequence of the anisotropic nature of magnetized plas- 
mas. This phase difference between circularly polarized modes 
results in a net rotation of the polarization angle by A$ /2. Fre- 
quently, this is extended to the case of elliptical plasma modes, 
which occurs when B is nearly orthogonal to the line of sight, 
or if the electrons are relativistic, producing Faraday conver- 
sion (also know as Faraday pulsations or generalized Faraday 
rotation), corresponding to oscillations between linear and cir- 
cular polarizations as well as a rotation of the polarization an- 
gle- 
Implicit in this description of the plasma processes under- 
lying Faraday rotation is the assumption that the two plasma 
eigenmodes do indeed propagate independently, or "adiabati- 
cally", i.e., there is no mode crossing. In the case of uniform 
plasmas (constant n e , B) this is true, and the foregoing expla- 
nation is complete. However, if the plasma is rapidly changing 
the situation is more complicated. For example, consider the 
propagation of the fast mode near the location of a magnetic 
field reversal, i.e., where k • B changes sign. Initially, the mode 
is nearly circularly polarized, with the electric vector rotating 
in the same sense as the gyrating electrons. As the compo- 
nent of the magnetic field decreases the mode becomes increas- 
ingly elliptical, changing into the linearly-polarized extraordi- 
nary mode when k • B = 0. On the other side of the magnetic 
reversal the mode again becomes elliptical and eventually cir- 
cular, though now with the sense of rotation reversed, following 
the sense of gyration of the electrons. Similarly, the slow mode 
reverses its polarization by evolving through the linearly polar- 
ized ordinary mode. The net Faraday rotation is then dependent 
upon J n e |B • d£\. Clearly the conditions for the applicability 
of Equation (1) depend upon adiabaticity being broken at field 
reversals. We call this the "standard" regime, with the regions 
where the adiabatic condition fails referred to as the "strongly 
coupled". However, the "super-adiabatic" regime, where the 
propagation is adiabatic throughout field reversals, can some- 
times be observed with interesting consequences. 

Since the failure of the adiabatic condition occurs near 
field reversals, we might expect that the standard and super- 
adiabatic regimes are distinguished by the properties of the 
eigenmodes when B • k = 0. Specifically, we might expect 
that when the eigenmodes change sufficiently slowly at field 
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FIG. 1 . — Mode behavior near magnetic field reversals. Top: The wave num- 
bers of the two plasma eigenmodes, with their quasi-longitudinal polarizations 
shown. Within the strongly coupled regime the amplitudes of the two modes 
effectively switch, producing the standard expression for the rotation measure: 
RMs oc J n e B ■ d£. Bottom: The terms in the adiabatic condition (Equation 2). 
Blue long dash, short dash and solid lines show (Ak/2) 2 , \Ak' /2\ and their 
sum, respectively. Magenta long dash, short dash and solid lines show x' 2 , <fi' 2 
and their sum, respectively. When the latter is larger than the former the prop- 
agation occurs in the strongly coupled regime, denoted by the shaded region, 
and significant mode conversion can occur. 



reversals, the propagation remains adiabatic. Near reversals 
the evolution of the eigenmodes are dominated by the in- 
crease in ellipticity, becoming significantly elliptical only when 
\6 — 7r/2| < Y = vg/v, where 8 is the angle between the mag- 
netic field and the line of sight and vb is the electron cy- 
clotron frequency. Thus, the eigenmodes change substantially 
on length scales on the order of Y(d6/d£)~ l near reversals. To 
ignore mode conversion we require that this be large in com- 
parison to Ak~ l , which at reversals is given approximately by 
c/nvXY 2 , where X = v^jv 1 in which v P is the electron plasma 
frequency 1 . Thus we arrive heuristically at the result that 
when cdBjdl <C TrvXY^ the propagation occurs in the super- 
adiabatic regime. 

More explicitly, as shown in the Appendix, the necessary 
condition for the local plasma modes to propagate adiabatically 
is 



Ak' 



> 



b' 2 + X ' 2 ; 



(2) 



where prime denotes differentiation with respect to £, <fi and 
X are the orientation and ellipticity angle of the polarization 
ellipse, respectively. However, to second order in the plasma 
parameters Ak is given by 



AA = 2tt—XY\ / cos 2 #+ (- 



sin 



(3) 



2 Equivalently, we require the energy of the mixing term associated with 
the evolution of the underlying eigenmodes be less than the energy splitting 
between the two modes 
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FIG. 2. — The growth the phase difference between the two plasma eigen- 
modes during propagation through a turbulent medium. Intervening cells of 
roughly uniform magnetic field with different orientations and strengths are 
shown on the top, with the source on the left and the observer on the right. The 
red line shows a path taken during standard Faraday rotation, with the phase 
difference accruing during solid segments and decrementing during dashed 
segments. The green line shows a similar path taken during super-adiabatic 
Faraday rotation, with the phase difference always accruing. Note the consid- 
erable enhancement of A<E> in the super-adiabatic regime, resulting in a corre- 
sponding enhancement in the rotation measure in this regime. 

where 9 is the angle between the magnetic field and the line of 
sight. When B _L k, Ak is reduced by a factor of Y /2 from its 
quasi-longitudinal value, and Ak' vanishes. On the other hand, 
near 9 = ir/2, the eigenmode ellipticity changes rapidly, giving: 



X 



Y sin6>(l+cos 2 6>) 



(4) 



4 (F/2) 2 + cos 2 6» ' 

where we have assumed Y <C 1 ■ If we further assume that the 
direction of rotation of the magnetic field is roughly isotropic, 
9' ~ 4>', implying that at reversals \' — <f>'/Y. Thus at the same 
time, the left side of Equation (2) falls substantially, the right 
side increases by a similar factor. This is illustrated explicitly in 
the bottom panel of Figure 1, which shows the left-side (blue) 
and right-side (magenta) of the adiabatic condition through a 
reversal. As a consequence, at magnetic field reversals the 
mode propagation almost always becomes strongly coupled, 
and in the limit of small Y the polarization propagates as in 
vacuum. 

A direct result of the failure of the adiabatic condition at 
magnetic field reversals is that the left-handed (for concrete- 
ness) fast mode maps onto the left-handed slow mode on the 
other side of the reversal, and vice versus. That is, there is an 
almost complete conversion of one mode to the other. Thus, 
we have the situation shown by the red line in Figure 2, as the 
polarized wave passes through regions of magnetic field with 
B • k > 0, A$ increases (solid lines), while when B • k < 0, A<£> 
decreases (dashed lines), yielding the well-known formula, 



A$ oc v 



n»B • d£ . 



(5) 



Note that the failure of the adiabatic condition at field reversals 
was critical to obtaining the standard result. 

Nevertheless, at sufficiently low frequencies, the adiabatic 
condition, Equation (2), can be satisfied at reversals. That is, 
because \Ak'\ vanishes and \' ^S> </>' where B • k = 0, in order to 
remain in the adiabatic mode propagation regime throughout, 
what we call the "super adiabatic" regime, we require 
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FIG. 3. — Mode behavior near magnetic field reversals in the super-adiabatic 
regime. Top: Same as in Figure 1, except now the modes are no longer cou- 
pled near B ■ k = 0. As a result, the fast/slow mode remains the fast/slow mode, 
though the local polarizations reverse, and RMsa °c J n e \h ■ d£\. Bottom: 
Same as in Figure 1, though note the absence of a strongly coupled regime. 



Note that up to a factor of two, this is precisely the condition we 
reached earlier by heuristic arguments. The critical frequency 
below which the propagation is super-adiabatic may be recast 
in terms of the local magnetic field reversal length scale, £g = 
ir/29' = 10 15 ^B,i5cm, and the local plasma parameters: 



(7) 



in which n and B are measured in cm" 3 and G, respectively. Al- 
ternatively, this may be written in terms of the total number of 
reversals along the line of sight, N ~ O'L/2-k, the rotation mea- 
sure determined at frequencies above i/$a, RMs = A~ 2 A$ ~ 
\~ 2 Ak\ g=Q L/ \/~N (where enters due to the cancellations dis- 
cussed previously) and the magnetic field strength: 



V B C 



2V27T 



1/2 



RMs 1 / 4 AT 1 / 8 



(8) 



13fi 1/2 RM s 1/4 AT 1/8 MHz, 



where in the final expression B and RMs are in G and radm -2 , 
respectively. Note that this is only a weak function of RM S 
and N, and thus primarily indicative of the local magnetic field 
strength in the source. As we shall see in Section 4, for a num- 
ber of sources this produces z^sa near 100MHz-l GHz. 

Since in the super-adiabatic regime there is no mode cross- 
ing at magnetic field reversals, the fast/slow mode maps onto 
the fast/slow mode, as illustrated in Figure 3. As a result, the 
phase difference between the modes only accrues, as seen by 
the green line in Figure 2. That is, for super-adiabatic Faraday 
rotation, 



A$ 



SA OC V 



\R-dl\ 



(9) 
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3. OBSERVATIONAL CONSEQUENCES OF 
SUPER- ADIABATIC FARADAY ROTATION 

There are three additional fundamental constraints upon ob- 
servable values of v>sa- The first is the requirement that within 
the Faraday screen, v$a must exceed the upper cutoff of the 
extraordinary mode: 



^sa > 



1 + W1+4- 



(10) 



Inserting the Equation (8) and assuming that there is consider- 
able separation between the cyclotron and plasma scale gives 
two conditions upon the magnetic field, 

B > 2.1 x 10- 5 RM s 1/4 A/- 1 / 8 4j 5 G, (11) 
corresponding to v SA > v P , and 

B <22RM S 1/2 AT 1/4 G, (12) 

corresponding to ^sa > v B 2 '■ Note that this produces lower and 
upper and bound upon 5 as a function of RMs. The former 
is shown for various values of N~ 1 I a Ib by the dashed lines in 
Figure 4. The latter is denoted by the upper-left greyed region 
in Figure 4. 

The second is that we must be able to ignore difference in the 
refraction of the two plasma eigenmodes. As the plasma eigen- 
modes propagate through the Faraday screen they necessarily 
follow different paths, dictated by their different indices of re- 
fraction. If the two paths pass through sufficiently distinct re- 
gions of the Faraday screen, Equation (1) is no longer valid, $ 
depending upon the differences in the plasma conditions along 
the paths in addition to the difference in the phase velocities of 
the two modes, and the RM is no longer constant. As shown in 
the Appendix, the condition that the modes propagate through 
strongly correlated plasma regions is v ^> vr where 



n l/3rl/6 

~ 61 , „ n l Wf' 3 kHz 



2.0D l 2 { 3 e B 2/ ^RM l/3 f l/3 MHz 



(13) 



in which / is the fractional variation in the plasma parame- 
ters, D = (Dq+Dj 1 ) - = D 2 2l0 22 cm where D and D s are 
the observer-screen and source-screen distances, respectively. 
The nature of this constraint depends upon the location of the 
Faraday screen. For example, when the rotation occurs in situ, 
D ~ L implying that D/£b — N, and giving vr that is a mod- 
erately weak function of N (oc N 1 / 2 or N 1 ^ 3 depending upon 
which of the above expressions one considers). For all of the 
sources we describe below, j^sa > vr, justifying our neglect of 
refraction for these objects. 

The third is the observational limits due to absorption in 
the ionosphere and heliopause (see, e.g., Bougeret 1996). 
The ionosphere limits ground observations to above roughly 
15 MHz, though down to 10 MHz is possible at excellent sites. 
A fundamental limit for space missions is the plasma frequency 
at the heliopause, estimated at roughly 9 kHz. These are shown 
by the hatched and solid grey regions, respectively, in the lower 
left. While a number of interesting objects lie within reach of 
ground-based radio astronomy, including nearby AGN and pul- 
sars, pushing to space-based, ultra-low frequency radio obser- 
vatories opens a new window upon a vast array of astrophysical 
environments. 

3 This also implies that we may ignore Razin suppression. 



If it can be observed, Faraday rotation in the super-adiabatic 
regime, what we term "super-adiabatic Faraday rotation", pro- 
vides a wealth of new information about the magnitude and 
line-of-sight geometry of magnetic fields. First, unlike stan- 
dard Faraday rotation the super-adiabatic rotation measure: 



RM SA = C / n e |B 



(14) 



does not suffer from a degeneracy between weak large-scale 
magnetic fields and strong tangled magnetic fields. Thus, 
RMs A is far better suited for the determination of n e and B via 
equipartition arguments than its high-frequency analog: 



RM S 



= C Jn e B- 



dt. 



(15) 



Second, generally, RMs A will be greater than RMs by an 
amount depending upon the number of magnetic field reversals 
along the line of sight. Typically, we may expect that looking 
through N independent turbulent cells, the excess phase differ- 
ence, and therefore RMs, to scale as A/" 1 / 2 , i.e., 



RM SA ~jV 1/2 RM s 



(16) 



Thus, measuring both RMs (at v > usa) and RMsa (at v < v$a) 
provides an estimate for the number of magnetic field reversals 
along the line of sight: 



RM 



SA 



RM S 



(17) 



However, this estimate may not be appropriate in the pres- 
ence of long-range order in the magnetic field, and more than 
a rough estimate requires a detailed analysis of the nature of 
the magnetic fields under consideration. For example, if there 
was some reason to not expect a randomly reversing field, but 
instead a band-structured field (e.g., as might be expected from 
pulsars or hydromagnetic turbulence within a given cell), then, 
RMs A /RMs ~ N. 

Third, as seen in Equation (8), the location of the super- 
adiabatic transition depends upon the local parameters of the 
plasma, though we have recast them in terms of the known 
RMs and B by assuming that the intervening Faraday screen 
isn't significantly clumpy. However, in many sources this as- 
sumption may not be appropriate. Measuring vs A provides a 
direct measurement of the local plasma properties at the mag- 
netic field reversals, and thus the distribution of the electron 
density and magnetic field strength along the line of sight. 
On the other hand, if we can assume that these quantities 
are smoothly distributed, measuring RMs A and j^sa provides 
a method to estimate the magnetic field strength, independent 
of the intervening electron density: 



B • 



(18) 



where i^sa and RMs are measured in MHz and radm" 2 , respec- 
tively. 

However, both bandwidth and beam depolarization are likely 
to be more severe problems for super-adiabatic Faraday rota- 
tion than its standard counterpart. Both are exacerbated sim- 
ply because RMsa > RMs. Beam depolarization, however, 
is also more likely due to the nature of super-adiabatic prop- 
agation. At field reversals the polarization is effectively re- 
flected about the magnetic field as the handedness of the po- 
larization eigenmode switches. The observed polarization is 
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therefore dependent in part upon the magnetic field direction 
in the last turbulent cell. This can produce large gradients in 
the polarization position angle on angular scales comparable 
to the turbulent scale. Furthermore, because probing the super- 
adiabatic regime requires pushing to longer wavelengths, it also 
generically results in larger beams. Nevertheless, for compact 
sources, those that subtend a small angular scale in comparison 
to the turbulence in the intervening Faraday screen, this effect 
is largely mitigated. 

4. IMPLICATIONS FOR SPECIFIC SOURCES: 
GROUND BASED OBSERVATIONS 

We now turn to the implications for specific classes of 
sources, and in particular estimate vsa- Figure 4 summarizes 
these. Here we dicuss those sources for which the super- 
adiabatic regime is accessible from ground based observations, 
i.e., u SA > 15MHz. 

4. 1. Ionospheric Faraday Rotation 

The ionosphere presents an opportunity to investigate the 
super-adiabatic regime via local radio experiments. The F- 
layer is roughly 300 km high with a typical free-electron den- 
sity of 10 6 cm" 3 . Using a typical terrestrial magnetic field 
strength of 0.5 mG, this results in an RM of roughly 4radm" 2 
at zenith, and potentially considerably more at oblique angles. 
This will be further enhanced at low frequencies by ionospheric 
refraction. Thus, while moderate in comparison to most astro- 
nomical sources, the atmospheric rotation measure is certainly 
detectable. 

Near the poles the Earth's magnetic field necessarily changes 
direction, with 6' ~ 7r/7?®. Thus, v^a — 40 MHz, which is well 
above the ionospheric cut-off. While ionospheric refraction is 
significant at this frequency, the differential refraction is suffi- 
ciently small (v R ~ 5 MHz). Thus, transpolar observations of 
the Jovian or Saturnian decametric radiation provides a natu- 
ral terrestrial radio experiment in which to study the transition 
from the standard to the super-adiabatic regimes in a well un- 
derstood plasma. 

4.2. Resolved Galactic Nuclei 
4.2.1. Sagittarius A* 

The supermassive black hole at the center of the Milky 
Way, Sagittarius A* (Sgr A*), is presently the best studied 
accreting black hole. Despite its comparatively small mass 
(4.7 x 1O 6 M , Ghez et al. 2008; Gillessen et al. 2009) and ex- 
traordinarily low luminosity (L/Ledd — 10~ 10 ), due to its prox- 
imity it has the unique distinction of being the only black 
hole to have been probed on sub-horizon scales (Doeleman 
et al. 2008). It also has the largest rotation measure observed, 
-5.6 ± 0.7 x 10 5 radm" 2 (Marrone et al. 2007; Macquart et al. 
2006). Estimates of Sgr A*'s magnetic field arise from both 
fitting simple accretion models to both the rotation measure 
results and the source spectrum (Agol 2000; Quataert & Gruzi- 
nov 2000; Falcke & Markoff 2000; Yuan et al. 2003; Marrone 
et al. 2007; Loeb & Waxman 2007), and range from 0.01- 
100G. The considerable uncertainty reflects the uncertainty in 
the location region responsible for the Faraday rotation as well 
as the geometry of the magnetic field along the line of sight. 
Nevertheless, these estimates imply that Sgr A* should enter 
the super-adiabatic regime between 30 MHz and 3 GHz. 

Typical length scales within the Faraday rotating medium 
surrounding Sgr A* are on the order of £ B ~ 10 3 GM/c 2 ~ 



10 15 cm. Even with the large RM's observed in Sgr A*, this 
still results in a v R ~ 1 MHz, and thus refractive decoherence is 
unimportant. 

4.2.2. M81* 

The supermassive black hole at the center of M81 is a 
second representative of a class of low-luminosity AGN. At 
7 x 10 6 M Q , it is moderately more massive that Sgr A*, and 
with L/Ledd — 10~ 5 , considerably brighter, though still consid- 
erably sub-Eddington (Brunthaler et al. 2006). The lack of po- 
larization in M81 at wavelengths above 3.6mm implies lower 
limits upon the value of RM of 10 4 radm" 2 and 4 x 10 5 radm" 2 , 
depending upon whether the source is beam or band depolar- 
ized (Brunthaler et al. 2006). Estimates for the magnetic field 
in the Faraday screen are again made by appealing to modeling 
the spectrum of M8 1 * and in analogy with Sgr A* (see above), 
giving a range of roughly 2-2 x 10 4 G, due to its larger lumi- 
nosity and mass. This corresponds to somewhat larger z/sa, 
near 3 GHz, though closer to the v sa — vb cutoff. Again refrac- 
tion is unimportant in this source. 

4.2.3. Other Resolved AGN 

Determining the magnetic field strengths within the Fara- 
day screen of nearby AGN other than Sgr A* and M81* is 
the primary difficulty in assessing the likelihood of observing 
super-adiabatic Faraday rotation in these sources. For other re- 
solved, sub-Eddington AGN we may estimate the magnitude of 
the magnetic field in the accretion flow via scaling arguments, 
equipartition and Sgr A*. Specifically, in terms of the luminos- 
ity in Eddington units, ^dd = L/L^ad oc M/MeAd, the electron 
(and therefore ion) density near the black hole scales as 

n e oc oc lEddM~ 2 M E dd oc 4ddM"' , (19) 
M z c 

where we have assumed that the radiative efficiency is not 
a strong function of mass for substantially sub-Eddington 
sources. The ion temperature is largely independent of black 
hole mass, reaching 10 n K in nearly all cases. Therefore, the 
equipartition magnetic field strength scales as 

B oc (n e kT) 1/2 oc 4dd M " 1/2 • (2°) 

We can also use this to estimate the expected scaling of rotation 
measures: 

RM oc n e Bdl oc 4dd M ~ 1/2 ■ ( 21 ) 

Note that these imply that brighter and lower-mass sources will 
have higher values of v% A . Since the rotation occurs locally in 

these sources, v$aIvvl oc 4dd^ 1//4 ' an d tnus refractive decoher- 
ence is unlikely to ever preclude observing the super-adiabatic 
regime for resolved AGN. 

The black hole at the center of M3 1 is shown as an example. 
Despite being much more massive, M ~ 1.4 x 10 8 M Q , it also 
is extraordinarily under-luminous, with £ E dd — 10" 10 , similar to 
Sgr A* (Li et al. 2009). The correspondingly rescaled RM and 
B is shown in Figure 4. More generally, polarization observa- 
tions of resolved nuclear regions will fill the upper-right corner 
of the B-RM plane, corresponding to typical z^sa ranging from 
the ground-based cutoff of 15 MHz to 10 GHz. 

4.3. Unresolved AGN 

Faraday rotation has been one of the primary methods for 
studying the magnetic geometry of AGN radio jets and cores. 
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FIG. 4. — Lines of constant fsA for various rotation measures and source magnetic field strengths, assuming N = 1 (though since fsA is so weakly dependent 
upon N it is unlikely to be significantly different in practice). We require that vja lie above the upper-cutoff frequency, which reduces to roughly the condition 
that i^sa > and > up (see Equations 11 & 12). The first of these corresponds to B < 22RMs '/ 2 Af~'/ 4 G, and is violated within the upper-left grayed region. 
The second reduces to B > 21RMs'/ 4 A' ~ l l*(}J 2 /xG and is shown by the thin dashed lines in the lower-right for £g = 10 15 cm, 10 18 cm and 10 21 cm. The lower-left 
grey-hatched region denotes where i/^a lies below the ionospheric cutoff (~ 15 GHz) and the low-left greyed region denotes where vsa lies below the heliospheric 
cutoff (~ 10kHz). The locations, and therefore the expected j^sa, f° r various classes of potential sources are shown, based upon observed RM's and estimated 
magnetic field strengths. Points and filled ellipses correspond to individual sources, while the empty regions are rough estimates of the envelope bounding the region 
each class of sources occupies. The represented source classes, and the sections discussing them in more detail, are: Sgr A* (4.2.1), M81 (4.2.2), M31 (4.2.3), AGN 
cores & jets and high-redshift radio galaxies (AGN & HzRG-i, 4.3), intrinsic rotation in Pulsars (Pulsars, 4.4), the solar corona (Sun, 5.1), the interstellar medium 
(ISM, 5.2), X-ray binaries (XB, 5.3), Galactic center region (GC, 5.4), Galactic molecular clouds (MC, 5.5), infrared/starburst galaxies (IG, 5.6) and the intracluster 
medium (Abell 1795 & HzRG-c, 5.7). 
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RM gradients have even produced direct evidence for a heli- 
cal field geometry in AGN jets (see, e.g. Kharb et al. 2009). 
Presently, rotation measures exist for more than 40 AGN, rang- 
ing from roughly 10 2 radm~ 2 to nearly 10 4 radm~ 2 (Zavala & 
Taylor 2004). These are usually determined via fitting VLB A 
observations at 8 GHz and 15 GHz, and in nearly all cases does 
not resolve the central source. In principle, these are associ- 
ated with the accreting black hole or jet, and thus are in many 
respects similar to the systems discussed in Section 4.2, though 
generally with considerably higher £ E dd- 

Estimating the magnetic field is again the primary challenge. 
We obtain an approximate magnetic field strength based upon 
an equipartition argument with the emitting nonthermal elec- 
trons, assumed to be accelerated in relativistic shocks either in 
the accretion flow or in the jet. In this case we expect the elec- 
tron power-law index to be roughly -2 between some minimum 
and maximum Lorentz factors, and L„ is given locally by 

L„~3 x 10 32 7 min n e B 3/2 y^- 1/2 ergs- 1 Hz- 1 , (22) 

where n e , B, V and v are given in cm" 3 , G, pc 3 and GHz, re- 
spectively (the flat spectrum of AGN are then produced by spa- 
tial structure in these quantities, see, e.g., Blandford & Konigl 
1979). The electron pressure is given by 

P e - \e e = ^'y min m e c 2 n e ln( ] (23) 

3 3 V 7min / 

If the ion pressure is similar, equipartition gives P e ~ (3B 2 /Stt, 
for typical /3's of 10-100. Thus, upon combining Equations 
(22) & (23), given v, L„, and V, we obtain a lower limit upon 
the field strength: 

B ~ 2 x 10 2 (A/py 2/1 Ll% ^/7y-2/7 mG , (24) 

where L v = 10 35 L„,35 ergs _1 Hz _1 , A = ln(7 max /7 m i n ) and the 
other quantities are given in the units described previously. 
Note that this is weakly dependent upon j3, and thus the sys- 
tem would need to be extraordinarily sub-equipartition for the 
magnetic field to be appreciably lower. For concreteness we 
choose A ~ 10 and (3 ~ 10 2 . 

For the AGN jets & cores for which we found RM's in 
the literature (Zavala & Taylor 2004; Kharb et al. 2009), at 
v = 15 GHz, L„,35 is typically of order unity. More difficult is 
estimating V, since the vast majority of AGN cores are unre- 
solved, and we can therefore only place an upper limit upon 
V (lower limit upon the equipartition value of B). To esti- 
mate this, for each AGN we inspected high-frequency VLBI 
images to determine approximate upper limits upon its angu- 
lar size, typically comparable to beam sizes of ~ 1 mas (Taylor 
1998; Lister et al. 1998; Taylor 2000; Lister & Smith 2000; Lis- 
ter 2001; Porcas & Rioja 2002; Zavala & Taylor 2003, 2004; 
Jorstad et al. 2005; Lister et al. 2009). Given the distances 
inferred from the source redshifts, this corresponds to limits 
upon the physical sizes of roughly 5 pc, or volumes of 10 2 pc 3 - 
10 3 pc 3 . The resulting estimates for the magnetic field strength 
are 3mG-30mG, and shown by the dark-blue region labeled 
"AGN" in Figure 4. 

In addition to the AGN discussed above, there are a small 
number of high-redshift radio galaxies (HzRG) exhibiting 
anomalously large rest-frame rotation measures. These in- 
clude OQ 172 (z = 3.535, RM~4x 10 4 radm -2 , Udomprasert 
et al. 1997), 3C 295 (z = 3.377, RM ~ 2 x 10 4 radm" 2 , Per- 
ley & Taylor 1991; Taylor & Perley 1992), SDSS 1624+3758 
(z = 0.461, RM ~ 1 , 8 x 10 4 radm" 2 , Frey et al. 2008), and PKS 
B0529-549 (z = 2.575, RM ~ 1.5 x 10 4 radm" 2 , Brodericket al. 



2007). These RM are an order of magnitude larger than those 
discussed previously, and thus may represent a distinct class. 
Their association with cooling-core clusters has led to specula- 
tion that the Faraday rotation is associated with the intracluster 
medium and not intrinsic to the broad-line region. Neverthe- 
less, if the rotation does occur within the core, this implies typ- 
ical magnetic fields of 0.5 mG, somewhat lower than those for 
other AGN due to their comparatively larger distances. This is 
shown in the dark-blue region labeled "HzRG-i". 

Resolved AGN exhibit strong RM gradients, frequently with 
RM dispersions that exceed the averaged RM. As a result, 
the RM for unresolved AGN cores is likely to be consider- 
ably smaller than the peak values. Furthermore, as we have 
stressed above, our estimates of the local magnetic field are 
highly uncertain. This is primarily a consequence of the unre- 
solved natures of most AGN, implying that our estimate of the 
equipartition magnetic field is only a lower limit. If the Fara- 
day rotation is occuring near the supermassive black hole, as 
appears to be the case in Sgr A*, then the AGN region should 
move up and to the right. Therefore, while for AGN typi- 
cal ^sa's presently lie just at the frequency cutoff for ground- 
based observations, in reality it may be closer to those for 
their resolved counterparts. Furthermore, for in situ rotation, 

v R ~ Q3N^ 2 e B l QRMl /3 MHz, still comfortably below v SA . 

On the other hand, if the Faraday screen is far from the region 
responsible for the emission, the intrinsic magnetic field may 
be closer to 1 ^G-IOO^G, typical of the interstellar medium. 
In this case, the i'sa's would lie well below the ground-based 
cutoff, though above the heliospheric cutoff, implying that in 
these sources super-adiabatic Faraday rotation may be observ- 
able from ultra-low frequency, space-based radio observato- 
ries. Thus, even just the detection of super-adiabatic Faraday 
rotation would determine the region responsible for the ob- 
served Faraday rotation. 

4.4. Pulsars: Intrinsic Rotation 

The rotation measures of pulsars can be quite large, reach- 
ing 3 x 10 3 radm" 2 in some sources. Traditionally, this has 
been interpreted to be due to Faraday rotation in the interstellar 
medium (ISM, see Section 5.2). Addressing the propagation 
of radio waves through pulsar magnetospheres has been stud- 
ied extensively elsewhere (see, e.g., Arons & Barnard 1986; 
Barnard & Arons 1986; Thompson et al. 1994), and is beyond 
the scope of this paper. Generally, it has been found that the 
propagation modes are quite different than those we have de- 
scribed here, due to structure within the magnetosphere and 
the presence of a dominant population of relativistic particles. 
Nevertheless, it is not clear that propagation through the ple- 
rion does not result in a substantial contribution to the observed 
RM's. 

Within the plerion, the wind is again dominated by relativis- 
tic particles, which have strongly elliptical polarization modes, 
complicating the determination of the expected RM's. In addi- 
tion, the wind has a relativistic bulk motion, which can signif- 
icantly alter Equation (1), (see, e.g., Broderick & Loeb 2009). 
For objects like the Crab, for which the bulk Lorentz factor can 
reach 10 6 , these relativistic effects can dominate the Faraday 
rotation, depending upon where the Faraday rotating medium 
is located. 

Given these caveats, observing super-adiabatic Faraday rota- 
tion in pulsars would present direct evidence for intrinsic Fara- 
day rotation in pulsars, with the attendant consequences for the 
geometry of pulsar magnetospheres and recalibrating Galac- 
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tic magnetic field studies. Given the uncertainties mentioned 
above, it is unclear what the relevant magnetic field strength 
is in this case. Here we take as a fiducial value that at the 
light cylinder, B ic ~ B smi (R NS /R lc f = (2nR m /Pcf, typically 
on the order of 10G-10 3 G, though in some cases consider- 
ably higher. While the magnetic field strength is a strong func- 
tion of radius (though a considerably weaker function in the 
wind region, where B is roughly proportional to r~ x , than in 
the magnetosphere in which B oc r~ 3 ), it does not reach inter- 
stellar values until roughly 3 x 10" 6 pc, vastly larger than the 
light cylinder. If the entire measured RM is due to the re- 
gion near the pulsar, the implied magnetic field strengths would 
place pulsars in the upper-left of Figure 4 for 749 objects, 
collected from Noutsos et al. (2008) and the ANTF catalog 
at www.atnf.csiro.au/research/pulsar/psrcat (Manchester et al. 
2005), shown in green and labeled "Pulsars". 

5. IMPLICATIONS FOR SPECIFIC SOURCES: 
SPACE BASED OBSERVATIONS 

Despite lying below the ionospheric cutoff, the super- 
adiabatic regime for many objects is potentially accessible 
from space (e.g., Jones et al. 2000). These necessarily have 
i/sa > 9 kHz, and are thus above the heliospheric cutoff. How- 
ever, in many cases they are considerably closer to their respec- 
tive refractive decoherence limits. Nevertheless, the observa- 
tion of the super-adiabatic regime provides a prime motivation 
to push towards a space-based, low-frequency radio emission. 

5.1. Solar Corona 

Faraday rotation measurements of polarized radio sources 
near the solar disk provide a unique method with which to 
probe the magnetic field of the solar corona. Within this re- 
gion the solar magnetic field is still quite strong, ranging from 
lOmG to 0.2G, and dense, with an electron density on the or- 
der of 10 3 cm" 3 . Nevertheless, the short path length results 
in only modest RM's, ranging between lOradm" 2 and zero 
(Mancuso & Spangler 2000). These imply that is SA ~ 3 MHz, 
shown by the yellow region labeled "Sun" in Figure 4. The 

refractive limit is roughly O^Djf^/fRM^MHz, and thus 
roughly an order of magnitude below this. Therefore, potential 
low-frequency, space-based radio observatories could directly 
probe the Corona for large-amplitude, long-wavelength Alfven 
waves, a leading contender for the mechanism responsible for 
heating the solar corona. Alternatively, observing the refrac- 
tive decoherence, characterized by a sustained departure from 
the A 2 -law, can place a direct limit upon the scale of the domi- 
nant wave modes. 

5.2. Pulsars: Interstellar Medium 

Traditionally, the observed RM of pulsars has been attributed 
to the propagation through the ISM. When coupled with a de- 
termination of the dispersion measure, Pulsar RM's provide a 
probe of the strength and distribution of the Galactic magnetic 
field (see, e.g., Men et al. 2008; Noutsos et al. 2008). Typical 
inferred magnetic field strengths for the ISM are on the order 
of l^G-10/iG, and place Pulsar polarization measurements 
in the orange region labeled "ISM" in the lower-left of Fig- 
ure 4 for 749 objects, collected from Noutsos et al. (2008) and 
the ANTF catalog at www.atnf.csiro.au/research/pulsar/psrcat 
(Manchester et al. 2005). For many pulsars, this extends 
above the heliospheric cutoff, reaching in some cases 100 kHz. 
However, less well defined is the refractive decoherence limit, 



which depends critically upon the turbulent scale of the ISM. 
That is, i/ R ~ lOD^^^RM-^kHz, and thus as long as the 
ISM varies primarily on scales above 30 pc refraction may be 
ignored. On the other hand, sub-parsec turbulence in the ISM 
could push f R above 200 kHz, rendering the super-adiabatic 
regime inaccessible. 

5.3. X-ray Binaries 

As the stellar mass analogs of AGN, X-ray binaries might 
be expected to produce sizable RM's. This is supported by our 
RM estimate for accreting black holes. With typical luminosi- 
ties and masses on the order of 10 _1 LEdd and 10M Q , respec- 
tively, Equation (21) implies tremendous rotation measures, 
RM ~ 8 x 10 21 radm" 2 . Therefore, Faraday rotation within 
the accretion flow will almost certainly be band-depolarized 
below the X-rays (i.e., for A < 3 A given a resolving power 
of 500). Despite this, the transition into the super-adiabatic 
regime doesn't occur until roughly 40 eV, making it all but im- 
possible to probe the accretion flow in this way. However, the 
magnetic fields and densities in X-ray binary jets are consider- 
ably lower, and there have been a handful of radio RM mea- 
surements of these features. 

At GHz frequencies, SS 433's jets are polarized at the 1% 
level, and have a rotation measure of roughly -119radm" 2 
(Gilmore & Seaquist 1980; Stirling et al. 2004). The mag- 
netic field strength in the region surrounding the jets, inferred 
from fitting the source structure at 5 GHz, is roughly 4 mG, con- 
siderably weaker than that associated with the core (> 20mG) 
though much stronger than the ISM field. 

During outbursts, GRS 1915+105 and XTE J1748-288 have 
both exhibited substantial polarizations, rising in both cases 
above 20% level. In GRS 1915+105, during these periods, the 
polarization angle of the ejected component evolved in a fash- 
ion consistent with rotation measures of 60radm" 2 -80radm~ 2 
(Miller-Jones et al. 2005). While it isn't clear what the strength 
of the magnetic field is in this source, estimates based upon 
equipartition arguments in the ejected jet give B ~ 1 mG. For 
XTE J 1748 -28 8, on the other hand, the rotation measure varied 
between 7radm" 2 -12radm" 2 , over a 10 day period (Brocksopp 
et al. 2007). During this time, the minimum energy required 
to produce the radio ejection was 2.6(c/kpc" 1 ) 8 / 7 x 10 42 erg, 
where the distance to XTE J1748-288, d, is less than 8kpc. 
With an volume of 9 x 10 48 cm 3 inferred from the outburst 
timescale, this implies an equipartition magnetic field on the 
order of 3mG, similar to those in GRS 1915+105 and SS 433. 

These three objects are shown in Figure 4 by the cyan re- 
gion labeled "XB". The combination of relatively low RM and 
mG fields again places j^sa — 1 MHz. Taking a typical scale of 
10 16 cm, this implies v R ~ 20D ^/gRM^ 3 kHz, and there- 
fore refraction is insignificant in XRB jets. 

5.4. Galactic Center 

The inner lOOpc of the Milky Way is an extraordinarily ac- 
tive region, exhibiting intense star formation (e.g., the Arches 
& Quintuplet clusters), numerous supernovae remnants and 
the mysterious non-thermal filaments (NTF's). Unique to the 
Galactic center, the NTF's are comparatively bright, highly 
polarized filamentary synchrotron sources extending in some 
cases for nearly 30pc. In addition to the large NTF's (the 
Southern & Northern Threads, the Radio Arc and the Snake), 
there are a number of shorter NTF's, extending for 5,pc-10pc, 
and in some cases misaligned with their larger brethren (Yusef- 
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Zadeh & Morris 1987; Lang et al. 1999; Law et al. 2008). 

The origin of the NTF's is presently unclear. The orienta- 
tion of the large NTF's appears to be aligned with a poloidal 
Galactic field, suggesting that they are associated with mag- 
netic flux tubes that carry the Galactic field through the central 
regions. This is supported by the large observed polarizations, 
up to 70% in some places, which require the presence large- 
scale ordered magnetic fields (see, e.g., Yusef-Zadeh et al. 
1997). However, this is in conflict with the extraordinary mag- 
netic fields required to maintain stability within the Galactic 
center's dynamic environment, nearly a mG (see, e.g., Yusef- 
Zadeh & Morris 1987). Furthermore, the presence of smaller, 
misaligned filaments appears to argue against a connection to 
large-scale Galactic fields (LaRosa et al. 2004). 

The RM for the NTF's varies considerably along the fil- 
ament, ranging from roughly 10 2 radm~ 2 to 6 x 10 3 radm" 2 
(Yusef-Zadeh & Morris 1987; Morris & Yusef-Zadeh 1989; 
Gray et al. 1995; Yusef-Zadeh et al. 1997; Lang et al. 1999). 
These are quite large in comparison to the RM's of roughly 
200radm" 2 (implying field strengths of roughly 20 fiG) ob- 
served for extragalactic sources that are viewed through the in- 
ner 170pc of the Galactic center (Roy et al. 2008), suggesting 
that the Faraday screen must lie close to, or even are intrinsic 
to, the NTF's. This is supported by the large RM gradients ob- 
served across some NTF's, which require turbulent lOmG ISM 
fields in the region, field strengths nearly an order of magni- 
tude larger than those inferred in the NTF's themselves (Yusef- 
Zadeh etal. 1997). 

As a result, similar to Pulsars and AGN, there is consider- 
able uncertainty regarding the region responsible for the ob- 
served Faraday rotation, and therefore the relevant magnetic 
field strength. Furthermore, given the dynamic nature of the 
Galactic center region it is additionally unclear if the field is 
ordered or turbulent. The red region labeled "GC" in the mid- 
dle of Figure 4 shows the best guess values for the NTF's the 
Northern Thread, the Snake, the Radio Arc and G 359.54+0.18, 
and roughly characterizes the uncertainty in the magnetic field. 
By comparison, taking a fiducial NTF transverse scale of 0. 1 pc 
and line-of-sight length of 3pc gives a typical v R ~ 40 kHz, 
well below j^sa- This is only weakly dependent upon the lo- 
cation of the Galactic center Faraday screen, rising to 0. 1 MHz 
for distances of up to lOOpc. 

5.5. Galactic Molecular Clouds 

Magnetic turbulence is thought to be critical for supporting 
structures within molecular clouds and generating the dense 
molecular cores that are the sites of star formation (McKee 
& Ostriker 2007). Despite this, there is still considerable un- 
certainty regarding the magnetic field strengths and geometry 
within individual objects, and therefore the role that magnetic 
fields play in the process of star formation. Typical values 
of the magnetic field in Molecular clouds range from Galac- 
tic value of a few /j,G to nearly a mG inferred from Zeeman 
splitting of molecular lines (see, e.g. Crutcher 1999; Han & 
Zhang 2007). The measured molecular lines are presumably 
produced in the densest regions, and therefore indicative of 
the fields nearest the sites of star formation. Nevertheless, the 
comparative lack of large magnetic fields in stars implies that 
at some point during the star formation process the magnetic 
fields must somehow dissipate (either by diffusion or reconnec- 
tion), and thus these may not represent the strongest magnetic 
fields present in molecular clouds. 

High-resolution maps of the rotation measures of molecular 



clouds have been obtained via 21cm observations of nearby 
star forming regions (see, e.g., Uyaniker et al. 2003; Wolleben 
& Reich 2004). These generally find considerable substructure 
and exhibit large RM gradients. The average average rotation 
measures are on the order of 30radm 2 , though with a disper- 
sion of roughly 3 x 10 2 radm 2 about the mean value, and reach- 
ing 10 3 radm" 2 in some places. This puts molecular clouds in 
the lower-left of Figure 4, denoted by the dark red region and 
labeled "MC", and i/ SA ~ 30kHz-2MHz, reflecting the uncer- 
tainty in the underlying physical conditions. 

As with other sources, the value of j^sa would all but de- 
termine the location of the Faraday rotating medium. Further- 
more, comparing RMsa and RMs would provide an indepen- 
dent measurement of the magnetohydrodynamic turbulence be- 
lieved to support these clouds. In particular, this would allow 
a comparison between the geometry of the magnetic field and 
that expected from the turbulent velocity field. 

As with other sources, uncertainty arises in v R associated 
with the uncertain scales of the location of and variations in the 
Faraday screen. Nevertheless, even for large molecular clouds 
(10 2 pc) and relatively small turbulent structures (0.1 pc) the re- 
fractive decoherence frequency is still only roughly 60kHz, al- 
most certainly well below the relevant estimates for fs A . 

5.6. Starburst Galaxies 

Infrared galaxies are presumed to be undergoing a phase of 
prodigious star formation. Therefore, one might expect that 
these offer the prospect of measuring the properties of ex- 
tragalactic star forming regions generally, and extragalactic 
molecular clouds specifically. Radio observations of infrared 
galaxies have indeed shown typical RM's of lOradm" 2 , though 
with a dispersion of roughly 50radm" 2 (though some reach 
nearly 200radm" 2 ), these are somewhat more narrowly dis- 
tribute about the mean (Heald et al. 2009). This is a reflection 
of the fact that we are necessarily averaging over a much larger 
region than we can resolve in Galactic molecular clouds. In- 
deed, the RM for the individual galaxies is comparable to that 
averaged over individual Galactic molecular clouds (Uyaniker 
et al. 2003; Wolleben & Reich 2004). 

Again there is considerable uncertainty in the strength of the 
magnetic field in the star forming regions of infrared galax- 
ies. We adopt a typical value of 10 yuG, though the dynamics of 
recent galactic interactions may raise this substantially. Never- 
theless, the region occupied in Figure 4 is between molecular 
clouds and the ISM, labeled "IG" and shown in magenta. This 
puts v^a — 10 2 kHz. This is, however, becoming uncomfort- 
ably close to the refractive decoherence scale, comparable to 
the 60kHz obtained in section 5.5. Thus, in this case even de- 
tection of the super-adiabatic regime will shed light upon the 
conditions within the star factories of starburst galaxies. 

5. 7. Intracluster Medium 

The magnetic field embedded in the hot, turbulent gas at the 
centers of galaxy clusters has an estimated strength of 10/iG- 
30yuG. Combined with the kpc path lengths, dense cluster 
cores are capable of produce substantial RM's. 

Indeed, the fact that HzRG's are located at the centers of 
cooling-core clusters suggests that the anomalously large RM's 
observed in these sources may be associated with the intra- 
cluster medium (ICM) and are not intrinsic to the galactic nu- 
clei. This produces estimates for the magnetic field 1-2 or- 
ders of magnitude lower than those obtained in Section 4.3, 
roughly on the order of 20 yuG. These are shown in the dark- 
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orange region labeled "HzRG-c", implying that in the rest- 
frame vsk — 1 MHz for these galaxies. Considering the high 
redshifts of these sources, this gives ^sa — 300 kHz. This is 
well above typical estimates of v R , which for a 20kpc core ex- 
hibiting large-scale structure on lOOpc scales is still roughly 
only 20kHz. 

Observations of typical radio galaxies embedded in the 
centers of clusters find rotation measures as high as 3 x 
10 3 radm" 2 (Ge & Owen 1993, 1994). Representative of these 
is J1346+268, located near the center of Abell 1795. Com- 
bined with electron density estimates from X-ray observations, 
this gives a lower-limit upon the cluster magnetic field strength 
of 20 /iG. This is roughly an order of magnitude lower than the 
ICM equipartition value, suggesting that the magnetic field is 
likely to be tangled on scales small in comparison to the clus- 
ter size. Producing the observed RM with a stochastic field re- 
quires roughly N = 10 2 , resulting in an order of magnitude dif- 
ference in RMsa and RMs. However, the corresponding ^sa 
is roughly 300 kHz, comparable to that found for the HzRG, 
assuming that the Faraday rotation occurs within the ICM in 
those sources. 

6. CONCLUSIONS 

Despite commonly being described in terms of the indepen- 
dent propagation of plasma modes, non-adiabatic effects at 
magnetic field reversals force the mode crossings that are criti- 
cal to understanding the standard expressions for Faraday rota- 
tion. For nearly all astrophysical systems, at frequencies above 
a few GHz, these effects dominate the propagation at magnetic 
field reversals. However, at sufficiently low frequencies it is 
possible to enter a "super-adiabatic" regime, in which the two 
plasma modes propagate independently throughout, resulting 
in rotation measures that are proportional to J n e \B\\ \dl instead 
of / n e Bud£. 

For large-scale ordered magnetic fields, there is little differ- 
ence between the standard and super-adiabatic regimes. How- 
ever, for magnetic field geometries exhibiting many reversals 
along the line of sight, the rotation measure in the super- 
adiabatic regime can be substantially larger than that deter- 
mined at high frequencies. Thus, comparisons between the 
standard and super-adiabatic rotation provides a unique means 
to measure the magnetic field geometry along the line of sight. 

The frequency at which Faraday rotation transitions between 
the standard and super-adiabatic regimes, v>sa, is itself depen- 
dent upon the local plasma properties at the magnetic field re- 
versals, though only weakly dependent upon the magnetic field 
geometry. Thus even the detection of super-adiabatic Faraday 
rotation, marked by a rapid increase in the rotation measure 
at low frequencies, provides localized information about the 
plasma density and magnetic fields. This has the potential to 
remove the degeneracy between propagation path length and 
plasma density or magnetic field strength, allowing an unam- 
biguous determination of the location of the Faraday screen in 
many systems. 

For known systems, vsa ranges from the heliospheric cutoff, 
approximately 10 kHz, to nearly 10 GHz, with regions contain- 
ing stronger magnetic fields and higher rotation measures tran- 
sitioning at higher frequencies. For resolved active galactic nu- 
clei, such as Sgr A*, M81* and M31, this lies above the iono- 
spheric cutoff at 15 GHz, and hence is observable from ground 
based radio telescopes. Some unresolved AGN, HzRG and Pul- 
sars may be observable from the ground as well. Detection of 
super-adiabatic Faraday rotation in any of these sources would 



provide striking confirmation of the presence of in situ Fara- 
day rotation. For X-ray binaries, infrared galaxies, molecular 
clouds, the Galactic center region, intracluster medium, inter- 
stellar medium and the solar corona, the super-adiabatic regime 
should be observable from space-based radio observatories. 

The search for super-adiabatic Faraday rotation is espe- 
cially timely given the proliferation of ultra-low frequency 
telescopes, such as the Murchison Widefield Array in Aus- 
tralia (80GHz-300GHz) and the Low Frequency Array in 
Europe (10MHz-240MHz), built ostensibly for the purpose 
of mapping 21cm emission from the high-z Universe, and 
the Frequency Agile Solar Radiotelescope in North Amer- 
ica (50MHz-21 GHz), currently under development (Bastian 
2003). All of these devices will collect the full compliment of 
Stokes parameters, and thus will be capable of searching for 
super-adiabatic Faraday rotation. In the far term, space-based 
radio observatories, such as Astronomical Low Frequency Ar- 
ray, have been discussed for more than two decades now (Jones 
et al. 2000). They provide the only means to access the 10kHz- 
15 MHz radio window, and beyond the confines of the Earth, 
can in principle do so with resolutions comparable to much 
shorter wavelengths at the Very Long Baseline Array . The abil- 
ity to directly probe the geometry of astrophysical magnetic 
fields along the line of sight should add to the list of science 
motivators for such a capability. 
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APPENDIX 
RADIATIVE TRANSFER REGIMES 
In general, Maxwell's equations give 

(V 2 - W + u?e) -E = 0, (Al) 

for an electric field E, and a dielectric tensor e. For plane waves 
propagating along the z-axis in a plane parallel medium (we 
address the possibility of refraction in the following appendix), 
this reduces to 

d 2 F 

_+u, 2 e-F = 0, (A2) 

dz L 

where F is the Jones vector (i.e., a two-dimensional vector 
constructed from the transverse components of E). For an 
anisotropic dielectric tensor, there will exist two nondegener- 
ate transverse modes defined such that 

Lu 2 eFi = kfFi. (A3) 

In the case of a plasma, these are in general elliptically polar- 
ized and aligned with the component of the background mag- 
netic field normal to the direction of propagation, i.e., 

F 1 =q( sin * ) and F 2 = q( C ° &X V (A4) 
^ V (cosx I V ~ lsm X J 

where the orientation of the polarization ellipses with respect 
to a set of axis fixed in space is determined by 

Q=( sin ^V (A5) 

I - sin <p cos J v ' 
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Then, F = F\F\ + F2F2, may be inserted into Equation (A2) to 
give 

F('+2is 2 <pF( + (k\ -p 2 -tp 2 + is 2 tp' + 2ic 2 tpip) Fi 

= 2(tp- ic2<£>) F 2 ' + (tp' - ic 2 tp' + 2is 2 (fip) F 2 

F 2 -2is 2 ipF 2 + (k\ -p 2 -tp 2 - is 2 (p' - 2ic2tptp) F 2 

= -2(tp + ic 2 ip) F[ - (tp' + ic 2 tp' - 2is 2 tptp) F\ , (A6) 
where a prime denotes differentiation with respect to z, c 2 = 
cos2x, 52 = sin2x, p = <fi' and tp = When tp = 0, these 
reproduce Forsterling's coupled equations (cf. Budden 1961; 
Ginzburg 1970). 

Thus far, no approximations have been made regarding the 
wave length or scale lengths of the plasma. From the form 
of Equations (A6), it is clear that if <p and tp vanish, the two 
modes will propagate completely independently. In the limit 
that tp and tp are small in comparison to k\ t 2, we may look for 
WKB solutions of the form 



F = 



J f k;dz 



(A7) 



and hence 



f[ + is 2 tpfi = (V- " ic2p)h c~^ Akdz 

f2-is2<pfi=-(ip + ic2'P)fie ilAkdz , (A8) 

where Ak = k\—k 2 terms on the order of tp 2 , tp 2 , tp', tp', ip<p, 
and were ignored as they are small by assumption relative 
to those that remain. Further expand the as 



/1 = u\6 



-if(Ak/2)dz 



and fi = U2& 



if(Ak/2)dz 



(A9) 



Then, 



Ak 
~2 
Ak 
2 



-S2P u\ = (tp-ic 2 p)u2 



u 2 + i 

which may be combined to give 

'Ak^ 2 
2 



-S2tp U2 = —{lp + iC2tp)U\ , 



(A10) 



♦1 



+ %l> 2 +tp 2 - 



s2tpAk- 



Akf 



— J +xp 2 +tp 2 -s 2 tpAk+i— 



u\ = 



w 2 = 0. (All) 



If the tp and tp terms are dominated by the Ak terms, then 

m,2 ^ const x^W, (A12) 
and thus the /; are constant. Therefore, in this limit the modes 
propagate independently (the so-called adiabatic regime). In 
the opposing limit, when tp and tp dominate Ak, then Equations 
(Al 1) are indistinguishable from the isotropic case (i.e., Ak = 
0), and therefore the polarization propagates unaltered (the so- 
called vacuum regime). This can be directly proved by solving 
for mi. 2 in this limit and expressing the answer in terms of F. 
Because the tpAk term will only be relevant when 

'AAA , , 

-\p~tp 2 + tp 2 , (A13) 

this term may be neglected in determining the limiting regimes, 
yielding 



tp 2 + p 2 < 



tp 2 + tp 2 > 



Ak 
~2~ 

Ak 



Ak! 



Ak! 



adiabatic 



vacuum. 



(A14) 



What remains is to identify these limits in terms of the phys- 
ical characteristics of the plasma. That v exceeds the upper 
cutoff of the extraordinary mode requires X + Y < 1, which in 
turn requires that the propagation occur above both the electron 
cyclotron and electron plasma frequencies. Under these cir- 
cumstances the ionic contribution to the dielectric tensor may 
be safely ignored. Thus, the k\ :2 are given by the standard 
Appleton-Hartree dispersion relations: 



fcl,2 = Wi/l- 



2X(l-X) 



2(l-X)-y 2 sin 2 6»Tr 



where 



T= yV 4 sin 4 6> + 4y 2 (l-X) 2 cos 2 <9. 
To second order in Y, this gives 



Ak^^-XYJcos 2 e+ f| 



(A15) 



(A16) 



(A17) 



The ellipticity angle is given by (see e.g. Ginzburg 1970; Bud- 
den 1964) 

(A18) 



tan x = x + sgn(x) Vl+x 2 



where 



As a direct result, 



Ysin 2 6 
2(1-X)cos0 



X = 



x 1 ^ Y (1 + cos 2 6») sin^ 
2(1 +x 2 ) ~ 4 (F/2) 2 + cos 2 6> ' 



(A19) 



(A20) 



where the finally expression assumes that Y -C 1 and X',Y' <C 
6'. Note that for ~ tt/2, x' - Y~ l 0' ^ F"V which at high 
frequencies (Y -C 1) is generally much larger than and <p' , and 
thus tp dominates tp near magnetic field reversals in the left- 
hand sides of eqs. (A14), and in particular is maximized at 
6 = tt/2. 

On the other hand, as long as X',Y' <C <f> near 9 = tt/2, 
I Ak' /2\ vanishes identically, due to the symmetry of AA:. This 
does not mean that |AA'/2| doesn't dominate (Ak/2) 2 else- 
where, which is clearly the case for the parameters shown in 
Figure 1. Nevertheless, sufficiently close to 9 = tt/2, the right- 
hand sides of eqs. (A14) are dominated by the {Ak/2) 2 term, 
and reach their minimum. 

Therefore, for the mode propagation to remain adiabatic for 
all magnetic field orientations (i.e., through an entire magnetic 
field reversal), it is both necessary and sufficient to have 



■0< 



Ak 
~2~ 



(A21) 



It is this condition that we use to define the "super-adiabatic" 
regime. Note that at magnetic reversals 



(A22) 



and therefore, much smaller than the dominant term in the 
above adiabatic condition, either tp 2 or (AA/2) 2 , as was pre- 
viously claimed. 



Ak 




Ak 


tp < max 






2 


0< 


2 
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BRODERICK & BLANDFORD 



Faraday 
Screen 




FIG. 5. — Diagram of the relevant angles and distances within the thin- 
screen approximation for the fast (green) and slow (red) modes. The maximum 
deviation is <5, located within the screen itself. 

REFRACTION AND FARADAY DECOHERENCE 

In the previous section we ignored the refraction of the radio 
wave. For propagation near the plasma and/or cyclotron fre- 
quencies this is not justified. Therefore, here we address the 
corrections and limits that refraction places upon the computa- 
tion of the Faraday rotation. 

Refraction potentially enters into the problem in two ways. 
First it provides a third scale over which the plasma eigen- 
modes change, though typically larger than 2ir/x' and 2ir/<fi', 
and thus not relevant. Second, and much more important, is 
the possibility that the two plasma eigenmodes refract differ- 
ently, propagating through significantly different regions and 
sampling distinct plasma conditions. It is this latter effect that 
we concern ourselves. 

Let us define Xi 2(77) and ki^iv) be the trajectories and wave- 
vectors of the two plasma eigenmodes, parametrized such that 
at rj = they begin at the source and at rj = 1 they arrive at the 
observer. The phase difference between the two modes is then 

/V* * fe .*»y (BD 

Jo V d7 l dr l J 

That this may be reduced to an integral along a single curve, 
e.g., x = (xi +X2)/2, requires that the two curves never separate 
further than the local plasma correlation length, roughly £g, 
leading to refractive decoherence. Note that in this case, it is 
not the total diffraction of the trajectories, but the differential 
trajectory that matters. Generally, this is a complicated global 
condition, requiring the complete construction of Xi 2 prior to 
addressing. However, in the "thin-screen" approximation, we 
can obtain a simple condition for when refractive decoherence 
may be safely ignored. 

Within the thin-screen approximation we assume that the 
Faraday screen's width is much smaller than both, the distance 
between the source and the screen, D s , and the distance be- 
tween the screen and the observer, Do- Within this approxi- 
mation, the trajectories are roughly straight lines, joined by a 



lensing angle at the screen. Geometrically, it is easy to see that 
Dscti = Doji and a, + 7, = /?,, from which it follows directly 
that, 

* B ir!^-<A-fl). (B2) 

L) s+ L) o 

That the maximum deviation between the paths taken by the 
two plasma eigenmodes is less than the characteristic plasma 
correlation length reduces to 6 <C £b, or 

IA-/3 2 |«^4 = |, (B3) 

L>sL>0 L> 

where D~ l = D~ s l +D3 1 . At this point, we must determine /?,. 

It suffices to consider propagation in the quasi-longitudinal 
regime since this dominates the differential refraction (because 
the difference in the indices of refraction is largest and the 
propagation nearly always occurs in this regime). Therefore, 
the dispersion relation is given by 

Di,2-\ [ki 2 -Lu 2 + uj 2 X(\±Y)\ , (B4) 

i.e., D\ 2 vanishes along the ray. The equations defining the ray 
are 

dX\,2 _ 0Di2 

di] 9ki 2 

(B5) 

dig ,2 _ dDi.2 _ uj 2 8X(1±Y) 
drj dx L2 2 dx\.2 

Outside of the Faraday screen, X = Y = 0, and k 2 = u 2 . In the 
weak deflection limit, we may treat the component of k, along 
the original direction as unchanged, the refraction simply in- 
troducing an orthogonal component. That is, we set 

/3 u2 ~-u 1 dri——x — ~--X(l±Y)f (B6) 

J dr) dxi, 2 2 

where / is the characteristic fractional variation in the plasma 
parameters on the scale of £ B . Therefore, for a Faraday screen 
of width Ib we obtain -/3 2 | — XYf. Should the screen 
be composed of may turbulent cells, the trajectories will dif- 
fuse, with \/3i-f3 2 \^ \/NXYf ~ y/L/e B XYf. Therefore, for 
the maximum deviation between the paths of the two plasma 
eigenmodes to be sufficiently small we require 

£ 3/2 

XYf<< W}T2- (B7) 

As with the super-adiabatic condition, we may view this phys- 
ically, as a lower limit upon £g, or observationally, as an upper 
limit upon v. In terms of the latter, we have 

(2 DL^V" 

the expression underlying those in Equation (13). 
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